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ABSTRACT 

We report the discovery of an ultra-luminous quasar J030642.51-I-I85315.8 (hereafter J0306-I-I853) 
at redshift 5.363, which hosts a super-massive black hole (SMBH) with Mbh = (1.07 ± 0.27) x 
10^° Mq. With an absolute magnitude M1450 = —28.92 and bolometric luminosity Lhd ^ 3.4 x 
J0306-I-I853 is one of the most luminous objects in the early Universe. It is not likely to be 
a beamed source based on its small flux variability, low radio loudness and normal broad emission 
lines. In addition, a z = 4.986 Damped Lya system (DLA) with [M/H] = —1.3 ±0.1, among the most 
metal rich DLAs at z > 5, is detected in the absorption spectrum of tbis quasar. This ultra-luminous 
quasar puts strong constraint on the bright-end of quasar luminosity function and massive-end of 
black hole mass function. It will provide a unique laboratory to the study of BH growth and the 
co-evolution between BH and host galaxy with multi-wavelength follow-up observations. The future 
high resolution spectra will give more insights to the DLA and other absorption systems along the 
line-of-sight of J0306±1853. 

Subject headings: galaxies: active — galaxies: high-redshift — quasars: emission lines — quasars: 
absorption lines — quasars: individual (SDSS J0306±1853) 


1. INTRODUCTION 

Supermassive black holes (SMBHs) reside in the cen¬ 
ters of most galaxies, with masses in the range of Mbh 
10® — 10^Mq. In the local universe, only a few objects 
have the measured black hole masses Mbh ^ 10^® 

(e.g. iMcConnell et al.ll2011l : iVan den Bosch et al.ll2012ll . 
Understanding when and how SMBHs form and grow 
is an important topic in extragalactic astronomy. One 
of the most direct ways to place constraints on the BH 
seeds and the growth of SMBHs is to study the masses 
of SMBHs hosted by luminous high redshift quasa rs (e.g. 
iVolonteri fc Reesll2006l : iTrakhtenbrot et aI1l2011ll . 

BH masses in quasars can be estimated by the 
single epoch virial mass estimate methods based on 
the low-ionization H/3 and Mgii emission lines (e.g . 
iVestersraard fc Peterson! 1^0061 : iMcLure fc Dunlonll2004J l. 
High- ionization lines such a s Civ can also be used 
(e.g. iVestergaard fc PetersonI 1200611 : however, CiV is 
thought to be less reliable than H (3 and Mgii in lim¬ 
ited S/N and resolution data fe.g. iSheii et al.l 120111 : 
Trakhtenbrot fc Netzerl[2012t iDennev et al.ll2013[l . 
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The Sloan Digital Sky survey fSDSS: lYork et al]l200n[l 
and the SDSS- Bl Baryonic Oscill ation Spectroscopic 
Survey (BOSS; iDawson et al.l I2013D have provided the 
largest high redshift quasar sample and resulted in the 
discovery of several hundred z: ± 4 qu asars with bolo¬ 
metr i c luminosity of s everal 10^® Lq (|Schneider et al.l 
120101: 1 Paris et ^ l2014|l . The SDSS quasar sample has 
been used to study the quasar luminosity function (QLE; 
[Richards et al.l l2006at IRoss et all 201311 and BH mass 

function IBHME: iKellv et al.ll201C : iKellv fc Shenl[20T^ 

at low and moderate redshifts. The BHME suggested 
that the maxi mum mass of a black hole in quasar to be 
^ 3 X 1O^®M0 (jKellv et al.ll201fll : IKellv fc Shenll2013[ l at 
those redshifts. From the combin yion of SDSS quasars 
and Stripe 82 faint quasar sample. iMcGreer et al.l (|2013[1 
presented the most detailed study on the z ^ 5 QLE, 
especially at the faint end. However, the slope of the 
QLE at the bright end still has large unce rtainties (e.g. 
[Richards et aLll^Obal : IMcGreer et ^11201311 . mainly due 
to the small number of quasars observed to have ex¬ 
tremely high luminosities (e.g. M1450 < —27.5). 

The discoveries of luminous z ^ 6 quasars (e.g . 
iFan et aT1[2001[ : rWillott et al.l[201(l[Bahados et akllToi^ 

and z > 6.5 qua sars le.g. iMortlock et al.l 120111 

iVenemans et al.ll2015ll with luminosity of ^ 10^^ Lq and 
MgII based BH mass of ~ 10® Mq, suggest that the fast 
growth of BH may start from very early epochs (z > 20) 
and involve massive se ed BHs {Mgeed > IO^Mq) (e.g. 
IVolonteri fc ReesI [2006fl . Recently, the first ten billion 
solar mass BH at z > 6 has been discovered in an ultra- 
luminous quasar J0100±2802 at z = 6.30, selected from 
SPSS and Wide-he ld Infrared Survey Explorer (WISE ; 
iWright et al.ll2010ll photometric data (Iwii et al.il2015[l . 
This discovery posts significant challenge to the theory 
of BH formation and the early BH-galaxy coevolution 
again. 
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Fig. 1.— The spectru m of J0306+1853. The red line shows the best fitting and the grey line shows SDSS composite spectrum 
(|Vanden Berk et al.ll2Q()ir) . The cyan line shows the noise spectrum. The blue points, from left to right, represent SDSS 2 -band, 2MASS 
J, H and Ks photometry. The blue arrow labels the position of a DLA at 2 = 4.986. The left inner plot shows the Voigt profile fitting of 

the DLA. Blue curve is the best-fit and gray shade includes 95% limit. The column density of the DLA is Nm = 10 ' - 0.12 cm~^. Two 

right small insets show the fitting of Civ and Mgii emission lines; the green lines represent power-law and Balmer continuum, the cyan 
lines denote Fell contribution and the blue lines are contributions from Civ and Mgii emissions. The gaps in the spectrum are due to the 
low sky transparency there. 


WISE mapped the wh ole sky at 3.4, 4 6, 12, and 22 
fim (Wl, W2, W3, W4). iWu et all (j2012fl suggested an 
efficient way to find high redshift quasars by combining 
SDSS and WISE photometric data, since most late-type 
stars have a significantly bluer W ISE colors than those 
of quasars. iCarnall et al.l (|2015ll used a similar method 
by combing optical and WISE colors to search z > 5.7 
quasars within VST ATLAS survey. We are conducting 
a ^ 5 quasar survey based on SDSS-WISE colors and 
have discovered several ultra-luminous quasars (E. Wang, 
et al. 2015, in preparation). Remarkably, we discovered 
an ultra-luminous quasar J0306-I-1853 at z = 5.363 with 
a bolometric luminosity of L^oi ^ 3.4 x and a 

BH mass of 1.07 x J0306-I-1853, along with 

JOlOO-l-2802, are the only two quasars yet known with 
Afi 45 o —29 and the BH mass up to ten billion solar 
masses at 2 ; > 5. Considering the high completeness of 
our selection method, these two quasars could be the only 
two ten billion solar mass BHs hosted by ultra-luminous 
quasars in entire the SDSS footprint. We also present the 
detection of a z ^ 5 metally-enriched Damped Lya sys¬ 
tem (DLA) in the spectrum of J0306-I-1853. Throughout 
this Letter, we use a A-dominated flat cosmology with 
Ho = 70 kms“^Mpc“^, Dm = 0.3 and Da = 0.7. The 
optical magnitudes in this letter are in AB system, and 
infrared magnitudes are in Vega system. 

2. OBSERVATIONS AND DATA 

J0306-I-1853 was selected as a high redshift quasar can¬ 
didate based on SDSS and WISE photometric data (F. 
Wang, et al. 2015, in preparation). J0306-I-1853 is un¬ 
detected in the u, g-bands but is relatively bright in 
the r, i, z bands (r = 19.89 ± 0.03, i = 17.96 ± 0.01, 
z = 17.49 ± 0.02) and is strongly detected in four 
WISE bands (Wl = 14.31 ± 0.03, W2 = 13.46 ± 0.04, 
W3 = 10.51 ± 0.10, W4 = 7.59 ± 0.16). With r -i = 
1.93, i — z = 0.47, J0306-I-1853 is in the area of the 
ariz color space str ongly contaminated by M dwarfs 
((Richards et al .1120 0211 . However, with Wl — W2 = 0.86, 


it is clearly separated from late-type stars. J0306-I-1853 
can also be selected by the Z-W2 /W1-W2 selection cri¬ 
teria used in ICarnall et al.l (1201511 . Besides, J0306-I-1853 
is also detected in Tw o Micron All Sky Survey (2MASS; 
iSkrutskie et al.l [200611 in J, H and Ks bands with J = 
16.50 ± 0.13, H = 15.70 ± 0.12 and Ks = 15.12 ± 0.13, 
respectively. The overall shape of the spectral energy dis¬ 
tribution (SED) of J0306-I-1853 built from SDSS, 2MASS 
and WISE photometry is consistent with t he type I 
quasar composite SED ([Richards et aLll^OObbll . 

The first low-resolution optical spectrum of this source 
was obtained with the Lijiang 2.4m telescope using 
the Yunnan Fainter Object Spectrograph and Camera 
(YFOSC) on November 25, 2013 (UT). A 22-minutes 
exposure was taken with the G3 grism (A/A A ~ 680 
at 7 OOOA) and a I'.'S slit. This low resolution spectrum 
clearly shows a sharp break at about 7700A, and is con¬ 
sistent with a quasar at redshift beyond 5.3. A possible 
z ~ 5 DLA system is also visible in the spectrum. To 
confirm this, a higher resolution optical spectroscopic ob- 
servation using Mage llan Echellette spectrometer (MagE; 
iMarshall et al. [200kl on the 6.5m Magellan/Clay Tele¬ 
scope, with 30 minutes integration time, was taken on 
January 3, 2014. Our configuration of the MagE opti¬ 
cal spectrum provides a resolution of A/AA = 4100 in 
the wavelength range from 6 OOOA to 10,000A with a I'.'O 
slit. The data were reduced by a custom-built pipeline 
(G. Becker, private communication) and flux calibrated 
using a standard star. 

In order to obtain the C iv and Mg ii based BH masses 
and the metallicity of the DLA system in J0306-I-1853, we 
took a near-IR s pectrum with the F olded-Port Infrared 
Echellette IFIRE: [Simcoe et aI1l2010ll . an IR echelle spec¬ 
trograph on the 6.5m Magellan/Baade Telescope. The 
observation was carried out on January 17, 2014, with an 
integration time of 60 minutes, and provided a resolution 
of i? ~ 5000 from 0.82 to 2.5 ^m. We red uced the data 
using the IDL based pipeline fFIREHOSE: lSimcoe et al.) 
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Fig. 2.— Distribution of quasar bolometric luminosities and 
black hole masses estimated from Mgii lines. The Red star rep¬ 
resents J0306-I-1853 and the light green c ycle represents the lu¬ 
minous quasar JOlOO-l-2802 at z = 6.30 dWu et a l.ll2015| h The 
brown points denote 2 : 4.8 quasars from iTTakhtenbrot et ahl 

(I2011J) . The blue squares denote SPSS high redshift quasars a t 
z 6 JKurk et al.l 120071 : IJiang et al.l 120071: IDe Rosa et al.l 1201111 . 
The gr een upper triangles denote CFHQS high -redshift quasars at 
z 6 dWillott et al.|[2010l : IDe Rosa et aI.II201iri . The purple lower 
triangle denotes the most dist ance quasar ULAS J1120-I-0641 at 
z=7.085 JMortlock et al.ll201in . Bl ack contours and g rey dots de¬ 
note SPSS low redshift quasars from IShen et al.l 11201111 (with broad 
absorption line quasars excluded). The error bars represent the Icr 
standard measurement errors, and the mean error bar for low red¬ 
shift quasars is presented in the bottom-right corner. The dashed 
lines denote the luminosity in different fraction of Eddington lu¬ 
minosity. Note that the systematic uncertainties (not included in 
error bars) of virial BH masses could be up to a factor of 0.5 dex. 

1201 in and the flux was calibrated with a nearby AOV 
telluric star. The magnitudes estimated from MagE and 
FIRE spectra are consistent with the SDSS and 2MASS 
photometric data (Am < 0 . 2 ) and thus imply a small 
variability of J0306+1853 in a long time baseline. Fi¬ 
nally, we combined the MagE optical spectrum and FIRE 
infrared spectrum and scaled it to match the SDSS z- 
band and 2MASS J, H and Ks-bands photometric data 
for the absolute flux calibration (See Figure 1.). The fi¬ 
nal spectrum yields a redshift of z = 5.363 ± 0.004 based 
on the Mg ii line fitting. 

J0306-I-18 53 is undetected in t he NRAO VLA Sky Sur¬ 
vey (NVSS: ICondon et al.lll99i^ . the completeness limit 
of which is about 2.5 mJy at 1.4 GHz. This constrains the 
radio lou dness of J0306-I-18 53 to be i? = ^ 76 

following iShen et all (120111) . where /ecm and /2500 are the 
flux density (f^) at rest-frame 6 cm and 2500 A, respec¬ 
tively. The radio loudness upper limit together with the 
small variability suggest J0306-I-1853 is not a beamed 
source. 

3. LUMINOSITY AND BLACK HOLE MASS 
MEASUREMENTS 

To derive line widths and continuum luminosities used 
for single-epoch virial BH mass estimators, we per¬ 
formed least-y^ global fitting to the combined optical 
and near-IR spectrum. The spect rum was de-re d dened 
for Galactic extinction using the iCardelli et al.l (|19890 
Milky Way reddenin g law and E(B-V) derived from the 
iSchleeel et al.l (|1998D dust map. The spectrum was then 


shifted to rest frame using the redshift based on Mg ii . 
We first fit a pseudo-continuum model to account for 
the power-law continuum, Fe ii emission, and Balmer 
continuum in the l ine free region. Templates of Fell 
emission were from iShen et al.l (|2nil[) by combining the 
iVestergaard fc WilkesI (|2001l) template (1000-2200A) and 
the lSalviander et all (l2007f l template (2200-3090A). The 
Fell template is broadened by convolving the template 
with Gaussians with a range of intrinsic line widths 
and then scaled to match the observed spectrum. We 
model the Balmer continuum following IDe Rosa et al.l 
(Unni), which assumed a partially optically thick gas 
clouds with uniform temperature Te = 15, OOOAT and nor¬ 
malized to the power-law continuum at Xrest = 3675A. 
From the best pseudo-continuum fit, we measured the 
rest-frame 3000 A and 1350 A power-law luminosi¬ 
ties as 2.04 X 10^^ ergs s“^ and 3.46 x 10"^^ ergs s“A 
and the absolute magnitude at rest-frame 1450 A as 
Mi 45 o,ab = —28.92 ± 0.04. By assuming an empir- 
ical conversion f actor from the luminosity at 1350 A 
(|Shen et al.l201lD . we estimate the bolometric luminosity 
as Lboi = 3.81 X Ti35q a = (1-32 ± 0.08) x 10"^® ergs s“^. 
We get the Eddington mass to be ^ 10 ^°Mq by assum¬ 
ing Eddington accretion, which is an appro ximate lower 
limit of the BH mass (e.g. iKurk et ^120071 ). 

We then fit individual emission lines after subtracting 
the pseudo-continuum. As two gaussian model provided 
a sufficiently good model of the line profiles of Mg ii and 
Civ given the S/N of our spectra, we use two Gaussians 
to fit Mg II and C iv emission lines, respectively. For 
other emission lines, we use one Gaussian to fit isolated 
fines and multiple Gaussian s to fit blended lines, same 
as did in IJiang et al.l (I2007D . The best fitting is shown 
in Figure 1. We noticed that a portion of red part of 
Mg II profile is missed in the FIRE spectrum due to the 
low transmission, which may cause some uncertainties in 
the estimation of Mg II line width. The FWHMs and 
equivalent widths (EWs) of Mgii and Civ are 5722 ± 
640 km s-\ 22.6 ± 2.8l and 7320 ± 680 km s-\ 25.1 ± 
2 . 1 A, respectively. 

To estimate the black hole ma ss of J0306-I-1853, we 
adopt the relation obtained by IVestergaard fc Osmeil 
(|2009D for Mg II estimator: 


Mb (Mg II) 

If)®-®® 

■ALa(3000 A)' 

0.5 

■FWHM(Mg II)' 

Mq 


ergs s“^ 


10® km 


and IVestergaard fc Peterso'r] (|2006l) for Civ estimator: 


Mb« (CIV) 6.66 

■ALa(1350 A)' 

0.53 

■FWHM(C IV)' 

Mq 

10 ^^ ergs s“i 


10 ® km s“i 

( 


These relations give the black hole mass of J0306-I-1853 
as MBH{Mg If) = (1.07 ± 0.27) x 10i° Mq and 
Mbh{C IV) = (2.16 ± 0.48) x lO^o Mq. The Edding¬ 
ton ratio is Ltoi/LEdd = 0.95 for Mgii based mea¬ 
surement and 0.47 for Civ based measurement. The 
uncertainties of BH mass is measured from the 68 % 
range of the distributions from fitti ng 200 mock spec - 
tra generated using the method in IShen et ^ (120111) . 
Note that it does not include systematic errors of us¬ 
ing different BH mass estimators. By using different 
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Fig. 3.— Velocity profiles of metal transitions for the 2 = 4.986 
DLA. The total column density of Si is Nsi = cm“^ and the 

metallicity of this DLA is [M/H] = —1.3 ib 0.1. See §4 for details. 

Fell te mplates (^e.g. lSalviander et al.ll2007tlTsnzuki et al. 
2006 ) and different estimators (e.g . IMcLure fc Dunlon 
2004: IVestergaard fc Osnieil l2009t iWang et al.l 120091 : 
Shen et al.ll2011l : |Trakhtenbrot fc Netzei]l2012ll . the Mg ii 
based BH mass of J0306+1853 could be in the range of 
7.3-20.7 billion solar mass and the systematic uncertain- 
ties of Mgll bas ed BH masses could be ~0.5 dex (e.g. 
iShen et 120131 1. However, the Mgii based BH mass 
measurement is still the most reliable method for mea¬ 
suring BH masses at high redshift. We will only use 
the Mg II based BH mass for subsequent discussion. Fig¬ 
ure 2 shows the distribution of quasar bolometric lumi¬ 
nosities and black hole masses estimated from the Mgii 
lines for a large sample of quasars from literatures over 
a wide r ange of redshifts. Adoptin g the Mbh estimate 
based on iMcLure &: Dunlod (120041 1 calibration, the BH 
in J0306-I-1853 is more mas sive by ~ 0.25 dex t h an th e 
most massive one found in iTrakhtenbrot et al 1 (i2(rrni 
or by 0.38 dex compared with the iDe Rosa et all (l2011f l 
study. It is close to the most massive bl ack holes at an y 
redshift, and together with JOlOO+2822 (iWu et al.ll2015ll . 
the most massive BH yet discovered at 2 ; > 5. Con¬ 
sidering the high completeness of our selection method, 
J0306-I-1853 and J0100+2802 could be the only two ten 
billions solar mass BHs hosted by ultra-luminous quasars 
at 2 ; > 5 in the SDSS footprint. 

4. A DAMPED Lya SYSTEM AT 2 ~ 5.0 

Damped Lya systems (DLAs) are atomic hydrogen gas 
clouds measured in absorptions to background quasars 
with a column density higher than 2 x 10^° cm“^, which 
are unique laboratories for understanding th e conversion 
of ne utral gas into stars at high redshift (iWolfe et al.l 
I 2 OO 5 II . However, the number of known high redshift 
DLAs are very rare and o nly about ten DLA s have 
been discovered at z > 4.7 (jRafelski et al.l[2?)T^ . The 
studies of metallicities of these high redshift DLAs sug- 
gest a rapid decline in m etallicity of DLAs at z ^ 5 
(IRafelski et ^IM^IMl . 

A z = 4.986 DLA system clearly presents in the 
absorption spectrum of J0306+1853 (Figure 1.). We 
determine the A^hi value of the DLA system by fit¬ 


ting a Voigt profile using the method of IRafelski et al.l 
(|2012ll . The best fit yields a column density of A^hi = 

lO^O-Soto.ia cm“^. We estimate the error on the A^hi 
value by manually selecting values that ensure that A^hi 
solutions with assuming a single component of DLA plus 
possible deblending Lya forest in the wings. However, 
we need to caution that the systematic effects of line sat¬ 
uration and continuum pla cement will lend larg er uncer¬ 
tainties in most cases (e.g. IRafelski et al.ll^l2l) and the 
systematic uncertainties for all the measurements could 
be ^ 0.1 dex. 

This DLA is also associated with a number of corre¬ 
sponding metal lines, including Cii A1334, Sin A1304, 
A1526, Oi A1302, Civ A1548, A 1 550 a nd MgiiA2796, 
A2803 (Figure 3). IRafelski et al.l (|2012ll suggests that 
S and Si are ideal elements to determine metallicity. 
We use the apparent o ptical depth method (AODM; 
iSavage fc Sembachlll991h to derive the total column den¬ 
sity of Si with Sill A1304, A1526 absorption lines, and 
find Nsi = cm“^. The rest-frame equivalent width 

of Sill A1304 and Sill A1526 are 0.070 ± O.OObA and 
0.10 ± O.OO 7 A, respectively. We estimate the metallicity 
of this DLA is [M/H] = —1.3 ± 0.1. Note that the error 
does not include systematic uncertainties. This makes it 
among the most metal-rich DLAs found to date at z > 5. 
Although the metallicity of the DLA in J0306-I-1853 is 
about 0.7 dex higher than the average metallicity of other 
z ~ 5 DLAs, the existence of such system is still con¬ 
sistent with_thfi_£a£iddecline in metallicity of DLAs at 
z > 5 (jRafelski et al.l[2?)l^ . 

5. DISCUSSION AND SUMMARY 

Although gravitational lensing is a possible explana¬ 
tion for its high luminosity and inferred BH mass, we do 
not expect a large lensing magnification as J0306+1853 
is an unresolved object in the FIRE acquisition image 
under the I'.'O seeing condition. Besides, J0306-I-1853 is 
not possible to be a beamed source as implicated from 
the small flux variability, small radio loudness and the 
normal broad emission lines. The bright-end slope of 
z ~ 5 QLF is measured to be /3 < —4, indicating a 
rapid decline in the space density of luminous quasars 
at high redshift (jMcGreer et al.ll20T^ . We searched the 
entire SDSS footprint (~ 14, 000 deg^) with our SDSS- 
WISE selection method and find J0306+1853 to be the 
only quasar with such high luminosity at this redshift. 
Since our selection method is highly complete for lumi¬ 
nous quasars at z > 5 (J. Yang, et al. 2015, in prepara¬ 
tion), J0306-I-1853 is possible to be the only quasar with 
Afi 45 o < —28.9 at redshift z > 5 in the SPSS footprint 
From the prediction using the z ^ 5 QLF (jMcGreer et al.l 
l2013f) . the number of quasars with such high luminosity 
4.7 < z < 5.4 is about two in the whole sky and the 
probability to find one such luminous quasar in the SDSS 
footprint with our selection method is ^35%. Together 
with J0100-(-2802 at higher redshift, J0306+1853 is the 
only other known quasar with M 1450 up to ~ —29 at 
z > 5. Although the discovery of J0306-(-1853 is con¬ 
sistent with the expectation from the current measure¬ 
ments of the QLF, its existence (and our finding of only 
a single object at this luminosity) shows that the lu¬ 
minosity function extends to much greater luminosities 
than previously explored, and sets a strong constraint 
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on the bright end slope. A detailed studies of the QLF 
will be presented in a forthcoming paper (J. Yang, et 
al. 2015, in preparation) based on a complete sample of 
S PSS—WISE se l ected z ~ 5 quasars. 

iKellv fc ShenI ()2013[ 1 estimated the BHMF and Ed¬ 
dington ratio function (BHERF) at z < 5 for Type 1 
quasars from uniformly selected ~ 58,000 SPSS quasars. 
Assuming there is no evolution from redshift of 4.75 to 
5.25, their BHMF predicts that there should be several 
tens quasars with BH mass large than 10 billion Mq at 
redshift between 5.0 and 5.5. Since J0306-I-1853 is the 
only quasar discovered in our survey with Mbh of 10^° 
Mq and Eddington ratio close to the order unity dis¬ 
covered, this would imply that there is a population of 
fainter quasars with lower Eddington ratio hosting lO^'^ 
Mq black holes. A complete survey of such objects re¬ 
quires high quality near-IR spectroscopy of a large sampe 
of fainter z ~ 5 quasars. 

Absorption spectra of high-redshift quasars provides 
one of the most powerful tools to study intergalactic 
medi um (IGM) in the early universe (e.g. iSimcoe et al.l 
1201111 . Future high signal-to-noise, high resolution spec¬ 
tra of ultra-luminous quasars as J0306+1853 will be very 
promising probes to the study of the IGM evolution and 
the shape of the ionizing background at z > 5. 
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